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Abstract. We study the volume of the nodal set of eigenfunctions of 
the Laplacian on the m-dimensional sphere. It is well known that the 
eigenspaces corresponding to E n = n(n + m — 1) are the spaces £ n of 
spherical harmonics of degree n, of dimension M ■ We use the multiplicity 
of the eigenvalues to endow £ n with the Gaussian probability measure 
and study the distribution of the m-dimensional volume of the nodal sets 
of a randomly chosen function. The expected volume is proportional to 
\/E n . One of our main results is bounding the variance of the volume 
to be O(^). 

In addition to the volume of the nodal set, we study its Leray mea- 
sure. We find that its expected value is n independent. We are able to 
determine that the asymptotic form of the variance is co J^ st . 
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1. Introduction 

Let M be a smooth compact manifold and / a real valued function on 
M. We define its nodal set to be the subset of M, where / vanishes and 
we are interested mainly in the nodal sets of eigenfunctions of the Laplacian 
on M. It is known [6j, that generically, the nodal sets are smooth subman- 
ifolds of M with codimension 1. For example, if M is a surface, the nodal 
sets are lines. One is interested in studying their volume (i.e. the length 
of the nodal line for the 2-dimensional case), and other nodal properties for 
highly excited eigenstates. It was conjectured by Yau that the volume of 
the nodal set is bounded from above and below by a multiple of the square 
root of the Laplace eigenvalue. The lower bound was proven by Bruning and 
Gromes [5] and Bruning [3] for the planar case. Donnelly-Fefferman's cele- 
brated result [7j resolved Yau's conjecture for real analytic metrics. However 
the general case of a smooth manifold is still open. 



1.1. Spherical Harmonics. In this paper, we study the nodal sets for the 
eigenfunctions of the Laplacian A on the m-dimensional unit sphere S m . It 
is well known that the eigenvalues E of the Laplace equation 

Af + Ef = 

on S m are all the numbers of the form 

(1) E = E n = n(n + m - 1), 

where n is an integer. Given a number E n , the corresponding eigenspace is 
the space £ n of the spherical harmonics of degree n. Its dimension is given 
by 



-n m ~ l . 



2n + m-lfn + m-l\ ' 

2 M = M n = ^ ~ ? , 

n + m — 1 \ m — 1 J (m— 1)1 

Given an integral number n, we fix an L 2 (S m ) orthonormal basis of £ n 

771(2;) = r/i(x), rj 2 (x) = t%(x), . . .,VAf(x) = rftr(x), 



giving an identification £ n = Mr 1 . For further reading on the spherical 
harmonics we refer the reader to [1], chapter 9. 
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1.2. Random model. We consider a random eigenfunction 



(3) f{x) = \J ^^^2a k r] k (x), 

k=l 

where a/% are Gaussian N(0, 1) i.i.d. which we assume to be defined on the 
same sample space f2. That is, we use the identification £ n = M.-^ to induce 
the Gaussian probability measure v on £ n as 

aV[I) € (2vr)^/2 ' 

where a = (aj) G R are as in Q. 

Note that t> is invariant with respect to the orthonormal basis for £ n . As 
usual, for any random variable X on f2, we denote its expectation EX. For 
example, with the normalization factor in ([3]), for every fixed point x G S m , 
one has 

irmi N 

(4) E[/(o;) 2 ] = ^^r ? ,(x) 2 = l, 

8=1 

a simple corollary from the addition theorem (see (|28p ). 
Any characteristic X(L) of the nodal set 

L = L f = {x G S m : f(x) = 0} 

is a random variable defined on the same sample space 0. We are inter- 
ested in the distribution of two different characteristics. The most natural 
characteristic of the nodal set Lf of / is, of course, its (m — l)-dimensional 
volume Z = 2(f). The study of the distribution of the random variable Z 
for a random / G £ n is one of the goals of the present paper. 
Berard [2] showed that the expected volume Ei? is 



EZ(f) = const ■ y/E n 
(see proposition II. 4p and Neuheisel [8] proved that as n — > oo, 

(5) V a r(Z) = (-^- T )=0 / L 



(m-l) 2 / I . - m ~ 1 

n 3m+l / \yV 3m+l 

Remark 1.1. Rather than taking a k standard Gaussian i.i.d., Neuheisel 
assumes that the vector a = (a^) G R^" is chosen uniformly on the unit 
sphere S^~ l . However, it is easy to see, that, since Z(f) = Z(b ■ /) for any 
constant b G R, both of those models are equivalent. 

The volume of the nodal line of a random eigenfunction on the torus 

T m = R m /Z m 

was studied by Rudnick and Wigman [11]. In this case, it is not difficult to 
see that the expectation is given by KZ(f Tm ) = const - y/~E. Moreover, they 
prove that as the eigenspace dimension J\f grows to infinity, the variance is 
bounded by 



VarZ(/ r ") = 



E 
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which, in particular, implies that the tails of the distribution of the normal- 
ized random variable ^ die. 

More generally, one may also consider a random model of eigenfunctions 
for a generic compact manifold M. Of course, for generic manifolds, one 
does not expect the Laplacian to have any multiplicities, so that we cannot 
introduce a Gaussian ensemble on the eigenspace. Let Ej be the eigenvalues 
and 4>j the corresponding eigenfunctions. It is well known that the Ej are 
discrete, Ej — ► oo and L 2 (M) = span{cpj}. 

In this case, rather than considering random eigenfunctions, one consid- 
ers random combinations of eigenfunctions with growing energy window of 
either type 

f L (x)= Yl a j<f>j(x) 
Eje[0,E] 

(called the long range), or 

f S ( x ) = aj4>j(x), 

(called the short range), as E — > oo. Berard [2] found that 

KZ(f L ) ~ c M ■ ^E 
and recently Zelditch [13] proved that 

EZ(f s ) ~ c M • VE, 

notably with the same constant cm for both the long and the short ranges. 

Berry [3] computed the expected length of nodal lines for isotropic, monochro- 
matic random waves in the plane, which are eigenfunctions of the Laplacian 
with eigenvalue E n . He found that the expected length (per unit area) is 
again of size approximately \JE n and argued that the variance should be of 
order log E n . 



1.3. Leray nodal measure. Another property of the nodal line we con- 
sider is its Leray measure (also called the microcanonical measure). Given a 
function / on S m , we define the Leray nodal measure to be 

(6) £(/) := limj-meas{x G S m : < e}, 

provided that the last limit exists. One may write the definition © of the 
Leray nodal measure formally as 

(7) C(f) := j 5(f(x))dx, 

S m 

where 5 is the Dirac delta function. 

As is well known, the limit Q exists when V/ ^ on the nodal set in 
which case 




dv'jx) 
:f(x)=o} |V/(x)|' 



where v' is the Riemannian hypersurface measure on the nodal set. This 
holds almost always on £ n (see section 12. ip . 
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The distribution of the Leray nodal measure on the sphere was also consid- 
ered by Neuheisel. As in case of the volume, one may compute the expected 
value 

EC ■ 



'2tt 

using a rather standard computation (see proposition II. 3p and Neuheisel 
proved that the variance is bounded by 

(8) V ar£ = o(^M = 0( ' 



Remark 1.2. Here, as in the case for the volume, Neuheisel considered 
a slightly different variation of the random model (see remark ll.ip . Even 
though the Leray nodal measure is not invariant under dilations, i.e. 

C(b-f) = - b C(f), 
those models are still equivalent asymptotically, as N —* oo. 

The Leray measure £(/' r " ) for the random eigenfunctions on the torus 
T n was considered by Oravecz, Rudnick and Wigman [10J. The expectation 
is given by 

EC(f T )- 



/2vr 

These authors were able to establish the variance to be asymptotic to 

Var£(/n~c.l 
for some c > 0, for m = 2 and m > 5. 

1.4. The expectation. 

Proposition 1.3. Forn sufficiently large, the expectation of the Leray nodal 
measure of the random eigenfunction is given by 

Proposition 1.4. One has 

(9) EZ(f) = c m ■ \[W n , 

with the constant c m defined by 

27T m/2 



(10) 



mT(f 



1.5. Statement of the main results. Our main results concern the vari- 
ance of the Leray nodal measure C and the volume Z of the nodal set. We 
improve on Neuheisel's results ([5]) and ([8]) , and need to use some of the steps 
in his work; however because some of the arguments in Neuheisel contain 
gaps, we need to redo them, partially accounting for the length of this paper. 

For C we were able to determine its asymptotics precisely. 
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Theorem 1.5. As n — > oo, the variance of the Leray nodal measure is 
asymptotic to 



, 2 m - 2 7r 21 f 2 r(f )|5 m | i 

< u > v " £ w~ — ( ro -i)! -jr 

One should compare the asymptotic result (jlip to Neuheisel's bound ([8]). 

Remark 1.6. Note that unlike the torus, our proof here works for any 
dimension m > 2, including m = 3, 4. The reason is that for the sphere, 
the so-called two point function u (to be defined, see (fT9|) ) is related to the 
ultraspherical polynomials, a standard family of orthogonal polynomials |12j . 
In particular, using Hilb's asymptotics for the ultraspherical polynomials, it 
is easy to show that the 4th moment of u is dominated by its second moment 
(see lemmas I A. II and lA.4|) . 

Unlike the spherical case, the two point function for the ci-dimensional 
torus is related to the distribution of points 

n = (m, . . . rid) £ 7L d 



so that 

||n||^ = nf + . . . + n 



•||2 2 .2 _ & 

4vr 2 ' 



For (i > 5 a strong equidistribution result for n implies in particular that 
the 4th moment of u is dominated by its second moment. For the two- 
dimensional case we used a special result due to Zygmund. The remaining 
cases d = 3, 4 are, to our best knowledge, open. 

Concerning the volume, we have the following result: 
Theorem 1.7. One has 

(12) VarZ(/) = O0| 

asymptotically as n — > oo. 

Note that theorem 11.71 implies that the variance of the normalized ran- 
dom variable Z := A with expectation 1, vanishes as n — > oo. Thus, in 
particular, the tails of the distribution of Z "die" , that is, for every e > 0, 
most of the mass of Z is concentrated in [1 — e, 1 + e] . In addition, theorem 
11.71 bounds the "typical" size of the tail of the distribution of Z (and thus 
of Z). One should compare (fT2l) to (pj|), obtained by Neuheisel. Partially 
motivated by the recent result [GW| for the analogous ensemble of random 
one dimensional trigonometric polynomials, it may be possible to improve 
the bound to E n /M . 

1.6. On the proofs of the main results. The spherical case offers some 
marked differences from that of the torus [10] and [11| . Unlike the torus, 
which is identified with the unit square with its sides pairwise glued, the 
sphere possesses a nontrivial geometry. In the course of the proof of the 
main results, one has to study the joint distribution of the gradients V/(x) 
and V/(y) as random vectors, where x,y £ S m are fixed, and / 6 £ n 
is randomly chosen. The main obstacle here is that for different points 
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x £ S m , the gradients live in different spaces, namely, the tangent spaces 
T x (S m ) which are, in general, different. 

One then has to canonically identify the spaces T x {S m ) via a family of 
isometries <p x smooth w.r.t. x. In reality, such a choice is not possible for 
every x, and we treat this complication in section l2~4l 

Once the geometric problems are resolved, the treatment of the so-called 
two-point (or, alternatively, covariance) function and its derivatives is more 
standard, related to the well-known ultraspherical polynomials (see |12] or 
appendix EJ). In particular, we find that the geometrical structure of the 
so-called singular sets on S m is less complicated than the singular set on the 
torus (see section 14. 2p . 

1.7. Some Conventions. Throughout the paper, the letters x, y and z will 
denote either points on the sphere S m or spherical variables and t will denote 
a real variable. For x, y £ S m , d(x,y) will stand for the spherical distance 
between x and y. The letters /j,, u, v will be reserved for measures, where 
the measure v will stand for the uniform measure on S m so that dv(x) = dx. 
Finally, given a set S, we denote its volume by \S\. For example, 

m+l 

(13) - 277 ^ 



' m+l - 



In this manuscript we will use the notations A <C B and A = 0(B) inter- 
changeably. 

1.8. Plan of the paper. This paper is organized as follows. Section [2] 
is devoted to the computation of the expected value of the Leray nodal 
measure and the volume, that is, proving propositions 11.31 and 11.41 where 
the rest of the paper focuses on the variance of those characteristics, i.e. 
proving theorems 11.51 and 11.71 The treatment of the variance in both cases 
will be divided into two steps. First, we express it in an integral form in 
section [3l We treat the integrals obtained in section [3] throughout section 01 
In case of the Leray measure, we will be able to give a precise asymptotic 
expression. In the case of volume, we give an upper bound. 

Appendix [A] will introduce the reader to the ultraspherical polynomials 
and will also provide all the necessary background we will need in this paper. 
The goal of appendix[B]is to prove that the set of "bad" (singular) eigenfunc- 
tions in the space of all the eigenfunctions, is "rare" in some strong sense. 
Finally, appendix O will prove a particular nondegeneracy result for the dis- 
tribution of the eigenfunctions and its gradients, needed to give meaning to 
the integral formula obtained for the variance of the volume given in section 
® 

1.9. Acknowledgements. The author wishes to thank Zeev Rudnick for 
initiating this research and for his help and support while conducting it. 
Many stimulating discussions with Mikhail Sodin, Dmitry Jakobson and 
Stephane Nonnenmacher are appreciated. The author is grateful to Sherwin 
Maslowe for proofreading this paper. I would like to thank CRM Analysis 
laboratory and its members for their support. Some part of this research was 
done during the author's visit to the Bielefeld University, supported by SFB 
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701: Spectral Structures and Topological Methods in Mathematics. Finally, 
I wish to thank the anonymous referee for his comments and suggestions. 

2. Expectation 

In this section we prove propositions 11.31 and 11.41 As a start, we wish to 
stay away from the set of the singular functions discussed in section 12.11 

2.1. The singular functions. In this section we define the notion of the 
singular functions and formulate the intuitive statement that they are "rare" . 
The proofs are given in appendix [Bl 

Definition 2.1. An eigenfunction / G £ n is singular if 3x G S m with 
fix) = and V fix) = 0. An eigenfunction / G £ n is nonsingular if V/ ^ 
on the nodal set. 

A nonsingular eigenfunction has no self-intersections. We denote Sing C 
£ n to be the set of all the singular eigenfunctions. First, we claim that as a 
set, Sing is "small". 

Lemma 2.2. The set Sing has codimension 1 in £ n . 

Now, given x G S m and b G R, we denote to be the set of all the 
eigenfunctions which attain the value b at the point x. That is, 

(14) VI = {fe£ n : fix) = b}. 

The set V§ is a hyperplane in £ n . 

Moreover, given G S m x S m and b = (&i, 62) S M 2 we denote 

(15) = {/€£*: /(*) = *!, fiy) = b 2 }. 

For x 7^ ±y, V^ ,y is an affine subspace of £ n of codimension 2, as it is easy 
to see from the addition theorem (see section |2~2|) . 

The next couple of lemmas establish the fact that the intersections of 
Sing with and Vu' y for x 7^ ±y, are of codimension 1. Lemma 12.41 is 
essential while treating the variance of the Leray nodal measure (section 

El). 

Lemma 2.3. For every x G S m and the set 

Singl := Sing n VI 

has codimension 1 in V^. 

Lemma 2.4. Ifx,y G S m and x 7^ ±y, then for every b = (61,62) G R 2 , fte 
sei Sing^' v := Sm^ n P^'^ /ias codimension 1 in 7^ 

The proofs of all the lemmas of this section are given in appendix [Bl 

2.2. Two-point function. We define the so called two-point function, also 
referred in the literature as the covariance function 

IOTII M 

(16) uix,y)= U %ix,y)=E[fix)fiy)] = L-i ^(^(y). 

k=l 

The addition theorem [jQ, page 456, theorem 9.6.3 implies that 

(17) u(x,y) = Q™icos d(x,y)), 
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where 

Q™:[-1,1]^K 

are the normalized ultraspherical polynomials defined and studied in appen- 
dix [A] Recall that d(x,y) is the spherical distance so that 

cosd(x,y) = (x,y), 

thinking of S m as being embedded into M m+1 . 

It is immediate that u is rotationally invariant, i.e. 

(18) u(Rx, Ry) = u(x,y), 

where R is any rotation on S m . In case y is not specified, it is taken to be 
the northern pole N G S m , that is 

(19) u(x) := u(x,N). 

For every t G [-1, 1], |Q™(i)| < 1 and \Q%(t)\ = 1, if and only if t = ±1. 
Therefore 

(20) (u(x,y) = ±l)^(x = ±y), 
and 

(21) (u(x) = ±1) 44> (x G {N, S}), 

where N and S are the northern and the southern poles respectively. 

2.3. Leray nodal measure. We will need the following definitions from |lUj . 
section 3. 

For e > 0, set 

£e(/) := — meas{x : |/(x)| < e} . 

so that £(/) = lim £ __ +0 >C £ (/). 
For a > 0, (3 > let 

£ n (a, 0) = {/ G £ n : |/(x)| < a => \Vf(x)\ > (3} . 

The sets £ n (a, (3) are open, and have the monotonicity property 

ot\ > a 2 =^ £ n (ai,/3) C £ n (a 2 ,(3) 

and 

Pi>02=>Sn(a,p 1 )C : £ n (a,P2) . 
Moreover, for any sequence a n , f3 n — > we have 

\ 5in^ = y £ n {a n , (3 n ) . 

n 

We have (cf. ITU], lemma 3.1) 
Lemma 2.5. For / G £ n ( a , fi) and < e < a, we have 

CeU) « 

where the constant involved in the ' <c' -notation depends only on a and (3. 
To prove lemma [231 we will need lemma 3.2 from [10] . 
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Lemma 2.6 (Lemma 3.2 from |10|). Let g(t) be a trigonometric polynomial 
on [0, 2tt] of degree at most M so that there are a > 0, (3 > such that 
\g'(t)\ > f3 whenever \g(t)\ < a. Then for all < e < a we have 

1 M 
— meas{t G [0, 2tt] : \g{t)\ < e} < — , 

where the constant in the ' <c' -notation may depend on m only. 

Proof of lemma UTR Let (0i, . . . m ) be the standard multi-dimensional spher- 
ical coordinates so that x G 5 m is parameterized by 

x = (cos 0i , sin 0i cos 02 > • • • > sin 0i . . . sin ip m ) 

for (0i, . . . TO ) G R := [0,7r] X . . . [0, ir] x [0, 2ir]. It is well-known that for 
4>i ^ 0,7r, 2ir, {^^} is an orthogonal basis of T x {S m ) and we have 





sin 0i • . . . • sin0 fc _i, 

£>9>jfe 

so that the Jacobian 

, , Dx 



J = J(0i 



; > fm; — 



31) • • • , <Pm; 

satisfies 

J = sin0™ _1 • sin0 2 m ~ 2 • . . . • sin0 m _i. 
Let < e < a. We write 

meas{x G 5 m : < e} 

as an integral 
(22) 



| J(0i, . . . , m )|(i01 



meas{x G S m : \f(x)\ < e} = J X (^)dx = J 

S m A,, 
in the spherical coordinates, where we denoted 

A e :={PeR: |/(P)|<e}. 
For P G R, 1 < k < m we define p k (P) = m I m a^ (-P). so that 

II S^fe II 

(23) ||V/(x)|| 2 =^+^ + ...+p^. 

We decompose 

A e = Wi u w 2 u . . . u w m 

with 

W fc := {P£A e : \p k (P)\ = max 
Note that on W k , 

\pt(p)\ > -L, 

by ([23]) and ||V/(a;)|| > j3 on A t . 
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Note that for 0^, k ^ ko fixed, g(4>k Q ) '■= /(</>i, ■ ■ ■ > 4>m) is a trigonometric 
polynomial in <j>k on either [0, ir] or [0, 2ir] of degree < n < y/E n with 
derivative 

d 



9'{4>) 



Pl(P) 



so that on Wk , |</(0)| > HaJ^H^/^;- Thus lemma |2~B1 implies 
meas{6* : \g(<f>k )\ < e } < n g "n • e. 

1 1 d<f> k() 1 1 

Therefore the contribution of VFfc to the integral (|22|) is 

| J|#l • . . . • # m < / meas{(/> : |3(^jfc )| < e}d</>i • . . . d(j> ko 



< f 1 1 1 1 • \/^n e #l • • • • #fto • • • #™ < eV^"- 



which concludes the proof of the lemma. 



□ 



We conclude the section with a formal derivation of proposition 11.31 A 
rigorous proof proceeds along the same lines as the proof of theorem 4.1 
in [10] (see section 4.2), using lemmas |2.5| 12.21 and 12.31 We omit it here. 

Formal derivation o f proposition ["Ol Given a function / € £ n , we write its 
Leray nodal measure formally as 



£(/) = J S(f(x))dx, 
s m 

see (J7]). 

Then, taking the expected value of both sides and changing the order of 
the expectation and the limit, we obtain 

(24) EC(f) = J E5(f(x))dx. 

Now, for each fixed x E S m , the random variable v = fix) is a linear 
combination of Gaussian random variables, and therefore, Gaussian itself. 
Its mean is zero and its variance is 1 by (J2|) . Writing the Gaussian probability 
density function explicitly, we have 

OO 1 2 

e-2 a , 1 



ES(f(x)) = E6(v) = 5{a)^=da 



— oo 

To finish the proof of this proposition we integrate the last equality on 
the sphere and substitute it into (|2"4"]) . 

□ 
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2.4. Choice of orthonormal bases for T x (S m ). For every x G S m we 
will need to identify 

(25) tf> x ■ T x {S m ) - R m , 

so that given a smooth function / on S m , the function 

V/ (X) G M m , 

is, under the identification (|25p . almost everywhere smooth (i.e, C fc , if / is 
C k+1 ) of argument x. 

Since we will be typically interested in the length of the gradient, we will 
require the identifications (|25p to be length preserving, namely, isometries. 
This is naturally accomplished, given a choice 

B x = {ef , . . . , e^} 

of an orthonormal basis of T x (S m ) for every x G S m , so that for every vector 
ei, all of its coordinates satisfy the appropriate smoothness condition. To 
do so, we consider the sphere without its southern pole S 

R:=S rn \{S}. 

Choosing an arbitrary orthonormal basis B = corresponding to the 
northern pole provides such a basis B x for every x G R by means of the 
parallel transport of B along the unique geodesic linking N and x on R. We 
choose an arbitrary orthonormal basis Bg of the tangent plane Ts(S m ) of 
S m at the southern pole. It doesn't affect any of the computations below, 
and we will neglect it from this point on. 

Let g(x) : S m — > R be any smooth function. We will use the notation 

d d 

de- 9{x) = d^ 9lx 
for the directional derivative of g(x) at x along ef , i.e. 

±g(x) = (Vg(x),ef). 

In case of ambiguity, i.e. if we deal with a two variable function 

h(x,y) : S m x S m — > R, 

we write -^xh or -^yh for the derivative of g as a function of x with y 

constant, or vice versa respectively. Similarly, we will use the notation 
V x g(x,y) G T x (S m ) and V y g(x,y) G T y (S m ) to denote the gradient of 
g(x, y) as a function of x or y respectively. 



Remark 2.7. Note that with the choice of the identifications (|25p as above, 
we have 

(26) V x d(x, y)\(x, N ) = -^yd{x, y)\(x )N ), 

which is going to be useful in simplifying the covariance matrix S (see section 

ED- 



Remark 2.8. In fact, for all our purposes, it is sufficient to make the choice 
of the orthonormal bases locally. Such a choice is possible for any manifold. 



DISTRIBUTION OF THE NODAL SETS 



13 



J (m+l)x(m+l) 



%rlm 



2.5. The covariance matrix, expectation. Given a point x G S m we 
consider the random vector (v,w) 6lx M m 

(v,w) = (f(x),Vf(x)), 

where we use the identification ([25]) . It is easy to see that being a linear 
transformation of a mean zero Gaussians, its distribution is a mean zero 
Gaussian as well. 

We claim that the covariance matrix of (v, w) is 
(27) 

Ef(x) 2 E[/(x)V/(x)] \ (\ 

v E[/(x)V/(x)] £ E[V/(x)*V/(x) 

First, 

(28) E/(x) 2 = it(cc,ar) = 1, 

by the definition flU]), ([El) of u(x, y) and (|20|). 
Next, we have 

E(/(a0V/0r)) = ^VE(/(x) 2 ) = Vl/2 = 

by @. 

Finally, we compute E[V/(a;)*V/(a;)] . For ? / j, we have 
def de 



E 



defde y - 
1 j 



u(x, y) 



0, 



x=y 



computing the second partial derivative explicitly in local coordinates (see 
section E3] for an explanation of the partial derivatives notations). 
For i = j, we have by the rotational symmetry on S m , 



E ( ( ^f /(x))2 ) = / E (w(*)-v/(*))«fc 



1 



m\S r ' 

E n 
m\S T - 



-E 



J Vf(x) ■ Vf(x)da 



1 



-E 



-E 



f(x) 2 dx 



E n 

m \S 



1 



m\S m 

E[f(x) 2 ]dx 



s r 

\2i 



f(x) ■ Af(x)dx 



m 



S m S r ' 
by the divergence theorem and ([!]). Thus 



(29) 



E(Vf(x) t Vf(x))=E(Vf(y) t Vf(y)) 



E„ 



m 



2.6. Riemannian volume. Let \ be the indicator function of the interval 
[—1,1]. For e > 0, we define the random variable 



2e 



X 



\Vf(x)\dx . 



Lemma 2.9 (Lemma 3.1 from [H]). Suppose that f G £ n is nonsingular. 
Then 

vol(r 1 (0)) = limZ £ (/). 
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Lemma 12.91 implies that the expectation of the volume and its second 
moments are given by the following. 

Corollary 2.10 (Corollary 3.4 from The first and second moments 

of the volume 2(f) of the nodal set of f are given by 

E(Z) = E(limZ £ ), E(Z 2 )=E( lim Z tl Z e2 ) . 

e^O ei,€2— >0 

Lemma 2.11. For every f 6 £ n and e > 0, one has 

Z e {f) = O(V^), 
where the constant involved in the 'O' notation depends only on m. 

To prove lemma [2.11|, we use lemma 3.3 from [11] . 

Lemma 2.12 (Lemma 3.3 from |11|). Let g(t) be a trigonometric polynomial 
on [0, 2tt] of degree at most M . Then for all e > we have 

1 

Ye 

{*:|fl(t)[<«} 



\g\t)\dt < 6M . 



Proof of lemma \2.11\ We write Z e in the multi-dimensional spherical coor- 
dinates (see the proof of lemma 12. 1 If) as 

7(01, 4>mY 



R 



V/(0 



1) • • • , ^mj 



|J|# 



1 • 



Note that in the spherical coordinates, for ^ 7^ 0, 7r, 27r 

1 0/ 



1, • • • , (PmJ 



k/ Kk<r 



in the orthonormal basis associated to { gf- } . Thus 



l|V/||-|J|«^ 



k=l 



Of 



Note that for cj)^, k 7^ ko fixed, f(4>i, ■ ■ ■ , (f> m ) is a 1-variable trigonometric 
polynomial in (j>k of degree < n <C y/E n . Therefore, 



fco=l 
by lemma [2.121 



3/(0i, . . . , m ) 



{ft[/(^l,...,0m)|<6} 



□ 



Now we are in a position to prove the main result of this section, namely 
proposition 11.41 

Proof of proposition \l-4\ We saw that 

E2(/)=Elim2 E (/) 
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by corollary 12.101 Lemma 12.111 and the dominated convergence theorem 
allow us to exchange the order of taking expectation and the limit to obtain 



(30) 



EZ(f) = limE2 £ (/). 



By Fubini's theorem, 
" 1 
Ye 



(31) EZ e (f) = E 



X 



fix) 



\Vf(x)\dx 



K^(x)dx, 



where K\{x) is defined by 
(32) K}(x):=E 



1 

Ye 



X 



fix) 



\Vf(x)\dv(f). 



We write K~ in terms of the random vector (v,w), introduced in section 
rfas 



w\\dfj,(v, w), 



where dfj,(v,w) is the joint probability density function of (v,w), namely 
mean zero Gaussian with covariance S given by (|27|) . Writing the Gaussian 
probability explicitly, we have 

1 



Kl(x) 



2e 



1 



X 



\w\\ exp 



1, \A dvdw 
-iv,w)E 1 (v,w) t ) = 

2 / (2vr)( m + 1 )/ 2 v / detS 



1 



— / exp ( — —v )dv I \\w\\exp 



1 ||?y|| 2 m\ m m l 2 dw 

2 En ) (2vr)(™+i)/ 2 £;r /2 



fc 

= YeI eXP (- \ V ^ dV ' ^ / M 6XP ( " ^' l|2 ) 

-e M m 

changing the variables 



dw' 



( 27r )(m+l)/2 ' 




n i 
W . 

m 



Following (f3"Tj) and ([3D]) , we integrate the last expression and take the 

limit e — ► to obtain 

EZ(f) = c m ^E~ n , 

where 



(33) 



m(2vr)( m + 1 )/ 2 



\w || exp 



1 



\dw 



Finally, substituting 



exp 



'l|2 



= V^(27r) 



m/2 _ 



rn+l 
2 



rlf 



^We use this opportunity to note that in [llj . page 7, the boundedness of 



^rX(~r^) is unnecessary, and, in fact, wrong. 
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(see e.g. [11J . page 7) and ([13]) into the last expression yields (fl~0j) . 

□ 

3. An integral formula for the second moment 

3.1. Covariance matrices, second moment. Similarly to the computa- 
tion of the expected volume, we will naturally encounter a random vector 
on 

Rxlx R m x R m , 

defined as 

(f(x)J(y),Vf(x),Vf(y)), 
for some fixed x, y G S m , where we again use the identification (I25p . We will 
use the rotational symmetry of the sphere to reduce the discussion to the 
case y = N is the northern pole. Thus we consider 

(34) Z := (v l ,v 2 ,w l ,w 2 ) = (/(*), f(N),Vf(x), V/(iV)) 
for some y G S m . 

It is obvious that the joint distribution of this vector, is mean zero Gauss- 
ian. It remains, therefore, to compute the covariance matrix. We need the 
following notations. 

Let D = D{x) be the vector in W 11 defined by 

D(x) = V x u{x,y)\ (x>N) G T x (S m ) - R m . 
Note that for x ^ ±iV, we may use fjlTf) to obtain 

D(x) =Q™'(d(x,N)) S m(d(x,N))V x d(x, y)\ (x>N) . 
It is then clear from (|26p that we then have 

(35) V x u(x, y) | {X:N) = -V y u(x, y) \ (a . jAr) . 
Finally, let 

H = H(x) = (hij) 
be the m x m matrix defined as 

(36) H = V x V y u{x, y) \ {X:N) , 
i.e. H = (hjk) with entries given by 

d 2 

kjk = de x -de y k U<yX,V ^ {x ' N) ' 

We will be in particular interested in the conditional distribution of 

Z 1 = (w 1 ,w 2 ) = (Vf(x),Vf(N)), 

conditioned upon f(x) = f(N) = 0. 

For the variance computation of the Leray nodal measure, we will need 
the distribution of the random vector 

Z :=(v 1 ,v 2 ) = (f(x)J(N)). 

It is distributed mean zero Gaussian as well. 

The covariance matrices of the random vectors above are given in the 
following lemma. 

Lemma 3.1. Let x G S. Then 
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(37) A 



1) The distribution of the random vector Z = (vx,v%) is mean zero 
Gaussian with covariance matrix given by 

1 u(x,NY 
u(x,N) 1 

(2) The covariance matrix of the random vector Z is the (2m + 2) x 
(2m + 2) matrix 

A B^ 
B l C / 

where A G M2x2 is given by (|3*7|) . B £ M2x2m is given by 

-D(x) 



B -\D (X ) u 
and C G M2 mX 2m given by 



(39) C 



— I m H 

m " L 

m " 



with the "pseudo-Hessian" matrix H = (hjk) ofu given by (I36p . The 
distribution of Z is nondegenerate for x ^ ±N (this is equivalent to 
S being invertible). 
(3) The covariance matrix of the conditional distribution of Z\, condi- 
tioned upon v i = V2 = is given by 



(40) n 



m 



H t J^nj i- u 2 \ —uD t D D l D 



We call the matrix £1 the "reduced covariance matrix" of Z\, and 
one has 

(41) detS = detvldetO = (1 - u 2 )detCX 

Proof. Part ([T]) of the lemma is evident from the definition of the two-point 
function. It is also clear that part ^ of the lemma implies part d3|), since 
one computes the covariance matrix f2 of the conditional distribution from 
S employing 



The nondegeneracy of the distribution of the random vector Z for x ^ ±y 
follows directly from appendixO The matrix T, is then invertible, being the 
covariance of a nonsingular joint Gaussian distribution. 

It remains, therefore, to prove part ([2]) of the lemma. It is clear that the 
block A is the same as the covariance matrix in part ([1]), i.e. given by ([3j 

Now by the definition, 

/E(/(x)V/(x)) E(f(x)Vf(N))\ 
[E(f(N)Vf(x)) E(f(N)Vf(N))J 

and we have already seen that 

E(f(x)Vf(x)) = 

in section 12.51 as well as 

E(f(N)Vf(N)) = 0. 
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Also 

E(/(iV)V/(x)) = V x E(f(x)f(N)) = V x u(x, y)\^ N) = D{x) 
and similarly 

E(/(x)V/(AT)) = -D(x), 
which finishes the proof of (|38|) . 



Finally, we compute C. By the definition 
C 



E(V/(x)*V/(x)) E(V/(x)*V/(JV))\ 
E(V/(AT)*V/(x)) E(Vf(NYVf(N))J ■ 

We have already computed that E(V/(x)*V/(x)) and E(V/(iV)*V/(iV)) 
are given by (|29D . Finally, 



E(V/(x)'V/(iV)) = V s V !/ E[/(x)/(iV)] = V^V^x, y)| ( , jiV) = ff, 
and similarly 

E(V/(iV)'V/(x)) = H*. 
This implies (I39D and finishes the proof of the lemma. 

3.2. Leray nodal measure. 



□ 



Proposition 3.2. The second moment of the Leray nodal measure is given 
by 

S m \ f dx 



U X 



|2' 



cm /■ 

(42) E£(/) 2 = L_J y 

where u(x) is the two-point function given by i\19h and (1170 . 

As in the case of expectation, we give a formal derivation of proposition 
13.21 omitting a rigorous treatment. A rigorous proof is obtained following 
the lines of the proof of theorem 5.1 in [lUJ (see section 5.3), using lemma l2~4l 
in our case. The convergence of the integral on the RHS of (|42jl . necessary 
to the proof, follows from (|50p and lemma B~T1 

Formal derivation of proposition \3.S[ We write the Leray measure as ([7]) 
again, so that 



(43) 



E£(f) 2 = E 



E[5(f(x))6(f(y))]dxdy 



S(f(x))5(f(y))dxdy 
S m S m S m xS m 

= \S m \ I E[S(f(x))6(f(N))]dx, 

S m 

by the rotational symmetry of the sphere. Now, for a fixed x E S m with 
x 7^ ±JV, the random variables v\ := f{x) and v% := f(N) are multivariate 
mean zero Gaussian with covariance matrix A given by f)3T[) . 
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Thus, writing the Gaussian measure explicitly, we obtain 
E(6(f(x))6(f(y))=E[6(v 1 )6(v 2 )] 

1 . , t da 



6(0,1)6(0,2) exp(— -aA a 



2 2vr\/deO 



2vrVdet A 2tt^1 - u(x,y) 2 ' 
Plugging this into (@2|) yields (|i2|) . 



□ 



3.3. Riemannian volume. 

Proposition 3.3. The second moment of 2(f) is given by 
(44) E(Z 2 ) = \S rn \ [ K(x)dx 



where 
(45) 

„, s 1 f exp(-±(wi,w 2 )n, 1 (wx,w 2 ) t ) dw 1 dw 2 

K(x) = — / lltuiim^l 



VI -u 2 JR^xR™ Vdetft (27r) w+1 ' 

where f2 = Sl(x) is defined by (I40p . 

Denote 

K eii€2 ( x ,y):=^-j £ ||v/(x)||||V/( 2 /)||x(^)x(^)^(/). 

To prove the proposition we will need the following lemma (cf. lemma 5.3 
in [IT]). 

Lemma 3.4. For (x, y) G 5 m x S m with x ^ y, one has the inequality 

E„. 



(46) K eue , 2 (x, y) « m - 

V 1 - 

where the implied constant depends only on the dimension m. 

The proof is almost identical to the proof of lemma 5.3 of [TT] 0. 
Proof. Write f(x) = (/, U(x)), where U(x) is the unit vector 



/ 1 Cm I 

U(x) = f 17 ±(r li (x)) i eS Af - 1 , 



where {rji(x)}^L 1 is the L 2 orthonormal basis of £ n chosen, and where we 
identify the function / with a vector in via ([3]). Note that 

(U(x),U(y))=u(x,y) 

is the cosine of the angle between U(x) and U(y). 
We have 

Vf(x) = DU-f 



^Note that in [11] . there is a misprint in the course of the proof of lemma 5.3. 
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where the derivative DU is a m x M matrix. Equivalently, 



(v/W), 



1 < i < to. 



By the triangle and Cauchy-Schwartz inequalities 

d 



||V/(z)||<£ 



due to 



by ©. 

Therefore 



c/(* 



9ei 



E 



5/ 



(47) 



K e e2 (x,y) < 



4eie 2 



_ E n 
m 

HI/1172 



I/(*)I<« 
l/(»)l<ea 



Consider the plane tt C R spanned by U(x) and U(y). The domain of 
the integration is all the vectors / £ M - ^ so that the projection of / on 7r 
falls into the parallelogram P defined by the perpendiculars 1% and to 
the endpoints of ±U(x) and ±U(y). Denote the angle a between the sides 
of P, computed as 

cosa = (U(x), U(y)) = u(x,y). 
We claim that the area of P is 

1 

area(P) = 4eie2 



V 1 - u ( x ^y) 2 

To see that, we assume, with no loss of generality that £2 cos a < e\ 
(otherwise exchange between x and y) and a £ (0, ^). Now if furthermore, 

£2 < ei cos a, 

then the line does not intersect the interval [0, eiU(y)}, and the sides of P 
are easily seen to have lengths and , and the angle between the sides 
of P is a, so that our claim follows. Otherwise (namely if £2 > eicosa), 
a little trigonometric computation shows that the lengths of the sides of P 
are again and -?^- and the angle between the sides of P is a. 

° sina sin a ° 

Write the multiple integral in (|47|) as the iterated integral 



(48) 



2 p -||/H 2 /2, 



df)dp, 



where the variable p runs over all the points of the parallelepiped P. The 
inner integral in (|48|) is 0(1). 

Indeed, note that for every f\ S 7T- 1 -, 



\\ p+ h f e -\lP+h\\V2 



2 + n/if: 



< (1 + 



e 
-ll/i 



pP+II/i 

2 /2 



2 )/2 
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since ||p|| 2 e U p H 2 / 2 is bounded. Our claim follows from convergence of the 
integral / k at-2(1 + ||iy|| 2 )e _ l' m ll l 2 dw. Therefore 

rWfW 2 / 2 df < area(P) « e x e 2 



l/(lf)l<«a 

Substituting the last estimate into ((171) proves (jl6|) . 



□ 



We give a formal derivation of proposition 13.31 Having lemma 13.41 in 
our hands, a rigorous proof of proposition 13.31 is identical to the proof of 
proposition 5.2 of [TT] and we omit it here. In the course of the proof one 
shows that 

K(x)= lim K eie2 (x,N) . 

ei,e2— >Q 



Therefore, taking the limit e±, e% — ► in (|46|) . we obtain 
Corollary 3.5. If u(x) 2 ^ 1 ffoera 

A-(X) « ; ^ = . 

Formal derivation of proposition Iff. 51 Corollary I2.1UI allows us to write an 
expression for the second moment formally as 



EZ{ff 



E 



5(f(x))\\Vf(x)\\5(f(y))\\Vf(y)\\dxdy 



S m xS 71 



and changing the order of taking the integration, we obtain 



EZ{ff 



E 



(49) 



S m xS r ' 



6(f(x))-\\Vf(x)\\-6(f(y))-\\Vf(y)\\ 



E 6(f(N))-\\Vf(N)\\-6(f(x))-\\Vf(x 



dxdy 
dx, 



by the rotational symmetry of the sphere. In fact, the integrand 



E 



5(f(x))-\\Vf(x)\\-5(f(y))-\\Vf(y)\\ 



depends on d(x, y) only (this is the isotropic property of the random ensem- 
ble E n ). 

Now for a fixed x £ S m with x ^ ±iV, the joint distribution of the random 
vector Z defined as in (|34f) is Gaussian with mean zero and covariance S = 
E(a>) as in lemma IBTTl Thus we may write 



E 



S(f(x))-\\Vf(x)\\-6(f(N)).\\Vf(N)\\^ 

5(vi) • \\wi\\ • S(vo) • WwoW exp( — (v, io)E -1 (v, wY) , 

V ' " " 1 ' " " V 2 V ' v ' ; (2^) m + 1 v / T^?v / detO^ 



substituting the explicit expression for the Gaussian measure and using ((IT 
(recall that O = fl(x) is defined by ((ID])). 
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Substituting v± = V2 = 0, we have 



E 



S(f(x))-\\Vf(x)\\-S(f(N)).\\Vf(N)\\ 
J \\w1\W\w2W exp(— -wT l ^ l w t ) 



dw 



(2vr) m +Vl -ttVdetn" 



To obtain the statement of the proposition, we integrate the last expression 
over S m and plug it into (Tl9j) . 

□ 



4. ASYMOTOTICS OF THE VARIANCE 

In this section we prove theorems 11.51 and 11.71 

4.1. Leray nodal measure. Here we use the ultraspherical or Gegenbauer 
polynomials (see appendix |A] for details). 

Concluding the proof of theorem \1.5[ Using proposition [221 (jl^P and propo- 
sition 11.31 we obtain 

(50) 

Var(Z) = ^ [ - dX ^ ^ [(- 1 <) ';- 

V 7 2tt J m ^/l-u(xY 2tt 2vr jf V^l-Q™(t) 2 / 

where fi = fj, m is the measure on / := [— 1, 1] defined by 

(51) = __.(! ^ 

It is easy to check that ^, = g*z^, where 5 : S m — > I is the function 

(52) sr(x) := cosd(x, N), 

and d is the spherical distance (recall that v is the uniform measure on S m ). 

Lemma 14.11 together with (I50p conclude the proof of the theorem once 
noting (J2J). 

□ 

Lemma 4.1. One has the following asymptotics 
l 



1 

1 



— 1 

where e(m; n) is given by 
(53) e(m; n) : 



d M (t) = 2 m - 2 vr m / 2 r(^)-i ZT + 0(e(m; n)), 



log n o 



n m , m > 3 
and /i is i/ie measure defined by (|5ip . 
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To prove lemma 14.11 we will divide the domain of the integral (i.e. the 
interval I := [—1, 1]) into two subintervals: B := [—1 + 1 — with cq 
constant, and B c := I \ B. We will show that the main contribution to 
the integral in (|50p comes from B, bounding the contribution of B c to that 
integral. 

We will reuse this partition while proving theorem 11.71 (see section |4~4"|) . 
This justifies devoting a separate section (namely, section [4^2]) to the treat- 
ment of B c . In analogy to the situation of [10] (cf. section 6.1) and |llj (cf. 
section 6.2), we will call B and B c the nonsingular and the singular intervals 
respectively. The proof of lemma 2J] will be finally given in section 14^31 

The singular and nonsingular intervals, as well as some of their properties 
will be given in section 14.21 The proof of lemma 14.11 will be finally given in 
section 14.31 

4.2. The singular interval. In the course of the proofs of theorems 11.51 
and 11.71 we are g°i n g to deal with the function 

1 



h(t) 



defined on [—1,1]. We wish to expand it into the Taylor polynomial of 
f[s) = s2 around s = as 

(54) = i + (o^Ql! + o((Q™(t)> 4 ). 

To be able to justify the expansion above, we will have to bound Q™(t) 
away from ±1, as in corollary IA.31 This corollary provides us with a subin- 
terval B C [—1,1] (which will be referred as the nonsingular interval) of 
large measure ^, such that Q™(t) is bounded away from ±1 for all t £ B. 
Giving a special treatment to its complement (referred as the singular inter- 
val, even though it is in fact a union of two disjoint intervals), we will show 
that its contribution is negligible (see sections 14.3.11 and I4.4.2j) . We give a 
rigorous treatment below. 

Let / be the interval I = [—1, 1]. Choose any < eo < 1 and the constant 
Co > guaranteed by corollary IA.31 corresponding to eo, assuming that n is 
large enough in the sense of corollary IA.31 We fix eo and cq throughout the 
rest of the paper and define the nonsingular interval 



B = B n := [ - 1 + \ 1 



CO -, COj 

n" n 2 



Corollary IA.31 implies that the expansion (|54f) holds on B with the constant 
involved in the 'O'-notation dependent only on eo- 
By an explicit computation, it is clear that 

(55) fi(B c ) « n~ m , 



where [i is the measure on / defined by (|5ip . 

Recall that /i is the measure on [—1, 1] induced from the uniform measure 
v on S m by g : S m — > [—1, 1] defined by (I52|) . We also define the spherical 
nonsingular set 

SB := g~ 1 (B), 
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and the spherical singular set 

SB C := S m \ SB. 

Since, as it was mentioned earlier, \i = g*v, it is evident that 

(56) u(SB c ) = fi(B c ) = 0{n- m ). 

The set SB is analogous to the nonsingular set in the sense of |10j (cf. 
section 6.1) and [TT] (cf. section 6.2). The structure of SB on the sphere 
(i.e., its projection B into [—1,1] by g) is by far simpler than that of the 
singular set on the torus, due to the lack of problems of arithmetic nature. 

4.3. Proof of lemma 14.11 

Proof. We write 

l 



-.dt 



+ 



L vi -mm u B 

This, together with lemmas 14.21 and 14.31 imply the result. 

4.3.1. The contribution of the singular interval B c . 
Lemma 4.2. One has 

dfx(t) 



□ 



(57) 



Vi - WW 



<C n 



Proof. We will bound the contribution of the integral on 

B c n[o, i] = [1-^,1], 

the rest being similar. Furthermore, we may assume by symmetry, that 
Qn (*) > so that 

1 1 



s/i^qxW Vi - QW) 

In what follows we will, consistently with appendix (Aj adapt the notation 

m — 2 



Writing t = cosip, we have (ft £ [0, ^] for some constant c\ > 0. Sub- 
stituting into Hilb's generalized asymptotic formula (see lemma lA.2p . we 
have 



Q™(cos^) = C 



tfi J Q (mft) 



sin^ (sin^) c 

for some constant C = C™, using the normalization defined by (|73p . Taking 
0, the value of the constant C is easily seen to be 



the limit 

(58) 

where 

(59) 



C 



lim 

4>->o 



J a {ncft) 



C = C r 



lim 



Jo 



<j>^0 (ft c 



n~ a C, 
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since Q™(1) 7^ (one can obtain an explicit expression for this constant 
using the expansion of the Bessel function into power series, see e.g. [9], 
page 57). 

Thus, the contribution of the singular interval to the integral, is, for n 
large enough 

(60) 

1 ci/n ci/ra 

(sin (A)" 1 " 1 r ^ m ~ 1 

< / , = =?dd> <C 



1-^ V V sin ^( sin *) a 



CI 

-in 



(1 + 0(^)2)^ + 0(^)2) 





ci 



l,m—l 



ib r ' 

= n- m / . ^ 

by ®. 

We claim that 

(61) l-C^l» Cl ^. 
Having (I61D proved would imply that 



'1— % Jo 

which is the statement of the lemma. 
To see ([6T]) , it is sufficient to show that 

1 — ^ ~ ;js 
lim — > 

i/>->0 -0 

and 

aw 



(62) x 

for every ip E (0, ci]. However the former inequality follows from the Bessel 
function expansion into power series around ■0 = (see [S], page 57, (9.09)) 

l-(7^ = a V> 2 + O(n 

for some constant ao > 0, so that the limit is positive. 

To see ([62j) . we note that in the course of establishing (f60|) . we showed 



Q-(cos^)-- JqW 
n 



Therefore, if (|62p is not satisfied, taking n large enough would contradict 

m(t)\ < 1. 

□ 
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4.3.2. The contribution of the nonsingular interval B. 
Lemma 4.3. 



(63) f 1 =-1 dfijt) =2 m ~ 2 7r m / 2 r(-)-— ^- + Q(e(m;n)), 

V ' J Jl - (Q™(tW 2 n™" 1 V V 



B 



where e(m;n) is given by (|53|) . 
Proof. On £? we may write 



Vi - (W)) 2 2 

(see section I4.2p . Integrating, we obtain 
1 



1 + Mp! + o(<OT) 4 ) 



= 5 / (Q™(t)) 2 ^) +o(J{Q™ 



(t)) 4 df,(t) 



B 

1 



VI " WW) 2 

1 

= 0{n{B c )) + (~ y (Q-(t)) 2 ^(t) + OCmCB 6 ))) + o(| {Q™(t)) 4 dn(t) 

-l -l 
= 1 (2 ^Wa r( ™)_^ + (n— )) + 0(e(m; n)) 

= 2m -2^2 r( |) . _i_ + 0(e(m . n)); 

as stated, by (f55|) and lemmas IA7T1 and E~4l □ 
4.4. Riemannian volume. The goal of this section is to prove theorem 

o 



4.4.1. Plan of the proof of theorem We have by proposition [37 

(64) Var(£(/)) = \S m \ ( K{x)dx - c m E n , 

JS m 

where 



,exp(— W2)Q 1 (wi,W2) t ) dw\dw2 



isT(a;) = — 7 === / ||«;i ||||w;2 



and c m is a constant given by (fTOj) . 

As in case of the Leray nodal measure, we divide the integration range into 
the nonsingular set SB and its complement SB C (see section [4~2]) , We bound 
the corresponding contributions to the integral separately (see lemmas 14.41 
and I4.5P . Using corollary 13. 5[ it is easy to relate the contribution of SB C 
to the last integral in (|50p . which we already bounded while treating the 
variance of the Leray nodal measure (lemma [4.2p . 

It then remains to bound the contribution of the integral on SB. Here we 
may write ^ 2 = 1 + 0(u 2 ) and one may show that, up to an admissible 
error, we may replace it by 1. We will define a new matrix S by 

n = ^(i-s), 

m 
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and notice that substituting S = into the integral, the identity matrix / 
recovers the square of the expected volume (KZ) 2 . Bounding the variance 
is then equivalent to "bounding" the matrix S in some average sense. 

To quantify the last statement we set o~(x) to be the spectral norm of the 
matrix S(x). We will show that the variance is bounded by 



To bound J o~(x), we use the trivial inequality cr(x) < Vtr S 2 . We will prove 
that J tr(S(x) 2 ) <C jj-, and together with the Cauchy-Schwartz inequality 
this implies the statement of the theorem. 

4.4.2. A bound for the contribution on the singular interval SB C . 
Lemma 4.4. One has 

K{x)dx <C E n e(m;n), 



SB C 

where e(m;n) is defined by (|53|) . 
Proof. We use corollary 13.51 to write 

f K{x)dx«. m E n f dx = E n f j=^f= <E n e(m;n), 

s { c s { c V 1 -^) 2 I V 1 -Qn(t) 

obtaining the last inequality by lemma FOl □ 

4.4.3. A bound for the contribution on the nonsingular interval SB. 
Lemma 4.5. 

1 2 E n 



J K(x)dx = -±-(E(Z)) 2 + 0( 



FN 

SB 

Proof. Define f^i = ilx(x) by 17 = — • The matrix Qi is symmetric, and 
positive for a set of x ^ ±N, since is such. Therefore it has a positive 
definite square root Pi = Oi. 

Intuitively, fii approximates the identity matrix I. To quantify this intu- 
itive statement, we introduce the matrix 
(65) 

m 1 ( D l D -uD l D - (1 - u 2 )^ 

b ~ 1 " E n 1 - u 2 {-uD'D - (1 - u 2 )^ D l D 

and its spectral norm a = o~(x), i.e 

a = max lad, 

l<j<2m 

where a^ are the eigenvalues of S. Note that, since fii is positive definite, 
S -C / in the sense that all its eigenvalues are in (— oo, 1). 
Changing the coordinates 

p 

w = \ — zPi, 
V m 
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we write the definition (|45p of K(x) as 

(66) K(x) = / IK^II-IK^allc-i"-" 1, <h 



my/1 — u 2 J (27r) 

l 2 ™ 



m+l 



where for a G M 2m we write (a)i € M m and (a) 2 G M m to denote either the 
first or the last m coordinates. 
We claim that 

(67) Pi = 1(1 + 0(a)). 



This follows from that fact that if S ~ diag(ai) then P\ ~ diag(\Jl — oq) so 
that 

Pi — I ~ diag(\Jl — oti — 1) <C diag(\a.i\) <C <rl, 



since \/l — y — 1 < \y\ on (— oo, 1). 

Moreover, by the definition of the spherical nonsingular set, on SB, u(x) 
is bounded away from 1, so that one may expand 

(68) l^ = i + ( u 2 )> 



,2 

where the constant involved in the 'O' notation is absolute 



J^+i / \\{z)i\\^\{z)2\\e-^\l + 0{u 2 )){l+0{a)?dz. 



Substituting (fBTl) and (f68l) into ([66]) . we obtain 

a 

m(2n) m+1 
Continuing, we have 

K ( x ) = -T^bTT I II Will • ||(^ 2 ||e-5(ll^ll 2 +ll^ll 2 )dz 1 dz 2 (l + 0( u 2 ) + 0(a) + 0(a 2 )) 



m(2iT 

slE, 



m+l 

in vlffm 



m(2TT) 2 



y Hz'He^ll 2 '!!^/) (1 + 0(u 2 ) + 0(ct) + 0(<t 2 )) 



* (EiJ) 2 (l + 0(u 2 )+0(<T) + 0(a 2 )), 



by © and (|33|) . 

Integrating on S"i?, we obtain 

K ( x ) dx ~ J^\( EZ ) 2 E n{^ + J u 2 dx + 0{v{SB c ))y 

SB SB 



by ([9]) and lemma |4TB"1 

To bound the last expression, we use (|56p . as well as, by the definition 
(|19p of the two-point function, we have 

l 

n 2 dx < J (Q™(cosd(x,N))) 2 dx = j (Q™(t)) 2 d^(t) « 1 

SB S" 1 -1 

by lemma [A. II and ([2j). This concludes the proof of the lemma. 

□ 

Lemma 4.6. For a fixed m, as n — ► oo, one has 
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(1) 



a(x) 2 dx — =. 

TV 



SB 

(2) 

Ja { x)dx« 
SB 



Jf 



Proof. Part [2] of the lemma clearly follows from part [T] by the Cauchy- 
Schwartz inequality. Thus we are only to prove part [H 

To prove the statement, we recall that a is by the definition the spectral 
norm of S, defined by (f65|) . To bound J a 2 , we use the trivial inequality 
a < tr(S 2 ). 

Since, by the definition of the nonsingular set SB, the two-point function 
u(x) is bounded away from 1, we may disregard the 1 — u 2 altogether. We 
define the matrix 

U)b ~{-uD t D-(l-u 2 )H t D t D 

We claim that 



J tr Sfdx = 0{n 5 - m ) 



This is sufficient for the statement of the present lemma, since then 

a 2 dx <C -=g / tlS idx <-^o { tr S 2 dx <C —r ■ n 5 ~ m <C -rj. 
h, n J h n J n M 

SB SB s m 

Now the elements of the matrix S 2 are bounded by elements either of the 
form 

du , Bu , Bu , Bu , 



the form 



or the form 



h h h U 



Bu Bu d 2 u 

\^N)- w \ {x ,N)-Q^r\ [x ,N) 



Be* 

11 ""12 ~~13~~H 



d 2 u , d 2 u 



Be? Be? ' 3p3 Be 2 

where in all the expressions above z may be either x or y (see section [2 
for an explanation of the partial derivatives notations). 

Using the Cauchy- Schwartz inequality again and the symmetry with re- 
spect to the variables, it suffices to prove the inequalities 

(69) J (||(*)) 4 ^«n 5 - 

and 

2 

,5— m 



f(—b 

J \BefBe2 



(70) 70^0*0 dx«n 
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We may compute the partial derivative in ([69]) (assuming x ^ ±N) as 

— Q™ (cos d(x,N)) = -Q™'( cos d(x,N)) sin d(x,N) — d(x,N), 
oe\ oe\ 

so that, since J^d(x, y) is obviously bounded on S m , it is sufficient to bound 

l 

(Q n 1, (cosd(x,N))smd(x,y)) 4 dx = J (Q%'(t)) A (l - t 2 ) 2 dfi(t), 

and thus (|69p follows from lemma IA.61 recalling the definition (|5ip of the 
measure \x. 

As for (|70|) . we write the second partial derivative in the integrand as 



-g^CooBdCx.Jv)) — 



sin d(x, N) — — d(x, N) 



so that, using a similar argumentation, we conclude that the integral in (|7Up 
is bounded by 

< J (Qn ''" '(cos d(x,N))f (sin d(x,N)fdx+ j (<2™'(cos d(x, N))) 2 dx 
1 1 

= / {Q r r(t)) 2 (i-t 2 ) 2 d^(t)+ J (Qr(t)fd^t). 

-1 -1 
Therefore, (]70p follows from lemmas IA.5I and lA.7i 

□ 

AAA. Concluding the proof of theorem 
Proof of theorem We write (I64p as, 

Var(Z) = j K{x)dx+(J K(x)dx-E(Z) 2 ^j < £ n e(m; n) + ^|= < Li ' 



IT 



SB C SB 

by lemmas 14.41 and 14, 5\ where we use 



E 

e(m;n) < -=, 



due to (P) and ©. □ 

Appendix A. Legendre and ultraspherical polynomials 

The ultraspherical (or Gegenbauer) polynomials P°(t) '■ [— 1,1] — > R of 
degree n generalize the Legendre polynomials P n (t) = P®(t). We use the 
corresponding normalized polynomials Q™(t) for an integral m > 0, which 
differ from P° (for a suitably chosen a), by a constant, defined by 

<r a) = i. 
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A.l. Definition and basic facts. The Legendre polynomials P n are the 
unique polynomials of degree n, orthogonal on [—1, 1] (w.r.t. the trivial 
weight function), normalized by P n (t) = 1. More generally, for a real 

a > —1, 

we define the ultraspherical polynomials -P"(i), being, up to a constants, 
the unique sequence polynomials of degree n, pairwise orthogonal w.r.t. the 
weight function on [—1, 1] defined by 

(71) u;(t) = (l-t 2 r. 

It is defined uniquely by the normalizing condition 

T(n + a + 1) 



(72) PZil) 



r(n + l)-r(a + l)' 



once we know that t = 1 is not a zero of P", see |12j . chapter 3.3. The 
ultraspherical polynomials is a particular case a = (3 of a more general 
class of polynomials, usually referred to as the Jacobi polynomials P% (see 
e.g. |12| for more information). 

While studying the spherical harmonics on the m-dimensional sphere, 
we are interested in the ultraspherical polynomials with a = ^j^-, and 
moreover, we would like to normalize it by setting its value at 1 to be 1. 
That is, we define 

P a (t) 

(73) (.)',: U > : " ' 



where 

(74) a := ™ 



For example, 



Ql(t) = P n (t) = P n °(i) 



are the usual Legendre polynomials. 

Throughout the section, we fix an integral number m > 2, and use the 
associated value of a, defined by ([74"]) . It is well known that is either 
even or odd, for the even and odd values of n respectively, and |Q™(t)| has 
a maximum at t = ±1. 

The function v = P"(t) satisfies the differential equation ( [12], page 60, 
(4.2.1)) 

(75) (1 — t 2 )v" — mtv' + n(n + m — l)v = 0. 

Due to its linear nature, it is also satisfied by v = Q™(t). The following 
recurrence relation ( [12], page 83, (4.7.27)) will prove itself as very useful 

(76) (1 - t 2 )P«'(t) + ntP^it) - (n + aji^t) = 0. 

Note that this recurrence relation is not satisfied by Q m (t) due to the dif- 
ferent normalization constants for P n and P n —\. 
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A. 2. Some basic results. Recall the definition (j5Tj) and (jTTj) of the mea- 
sure fi = fi m and the weight function u m respectively. We note that 

27T m / 2 

so for purposes of giving an upper bound only, we may disregard the differ- 
ence between d[i and ujdt. 

Concerning the 2nd and the 4th moments of the ultraspherical polynomi- 
als, we have the following: 

Lemma A.l. For m fixed, the second moment of the normalized ultras- 
pherical polynomials is 
l 

<r(f) 2 <i„(*) = 2"" V' 2 r(?) ■ -Lj + 0(, 



J - II'- 

-1 



as n — > oo. 

Proof. One has ( [12J, (4.3.3)) 

(77) / p:(tMt)dt = — — \ n , 2 > — -, 

y ' J n\J \J 2n + m- ir(n + ir(n + m-l 

-l 

and (|72p implies that 

(78) P m (l) = { 2> c-n a 

{ ' a { ' r(n + i)r(f) 

Thus, using the definition ([73]) of the normalized ultraspherical polynomials, 
we obtain 

2 2™' 1 r(n + l)r(f ) 2 2vr m / 2 



Q^(t)^ M (t) 



2n + m- 1 T(n + m- 1) T(f) 
2 m 7r m / 2 r(f) n! 



2n + m — 1 (n + m — 2) ! 

= r irt(5 • i • — *-^(l + O(-)) 
v 2 ; n n m " 2 ^ V y 

= 2 m - 1 7r m / 2 r(^)^i T + 0(n- m ), 
2 n m 1 

as stated. □ 
Lemma A. 2 (Hilb Asymptotics (formula (8.21.17) on page 197 of Szego)). 

1 /2 

(79) (1 sing)»P»(co S 0) = jV-" r(n + " + 1} ( ^ J a (N0) + 6(O), 
2 n! ysmfc' y 

uniformly for < 9 < vr/2, where N = n + J a is the Bessel J function 
of order a and the error term is 



5(9) < 



jV/2 (n- 3 / 2 ), en' 1 <9<tt/2 
[9 a+2 0(n a ), < 9 < cn- 1 . 
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Remark: It is clear, that 

(80) n _ a r(n + a + l) =1+ 1 

n\ n 
so we will usually omit this factor. 

Corollary A. 3. For every eo > 0, there exists a constant cq > depending 
on m only, such that if c > Co and t € [0, 1 — %], where n is large enough 
so that the interval above is not empty, one has 

\Qn(t)\ < eo- 

Proof. Let t = cosd. Then if < t < 1 - 9 > C ■ ^ for some absolute 
constant Co > 0. Lemma IA.2I implies that one has 

\P:(cos6)\ <C-±-\J a (N9)\ 

olll U 

for some absolute constant C > 0. We bound it by 

n a 

|P«(cos0)|<C„|J a (iV0)|, 

c 

so that (|78p implies that 

\Q™(cos9)\<-^\J a (N9)\<e , 
c o 

provided that we choose Co large enough, since J a is bounded. □ 
Lemma A. 4. The 4th moment of the ultraspherical polynomials satisfies 



Q™(i) 4 d/z(t)«e(m;n), 



where e(m;n) is defined by (|53[) . 

Proof. We will limit ourselves to the interval [0, 1]. To prove the statement 
there, we invoke the generalized Hilb's asymptotics (lemma |A.2|) , 
We have, using ([80]) . that 



(sin 9) m - 2 (P: (cos 9)) 4 « 4 (iY0) _^_ + _i!g)_ 



(sin6») T 



sm 



and claim that 

i 

(81) / P%(t) 4 dti(t) « 



— Sj-, m = I 
ra m " 4 , m > 3 



-i 



We have 



tt/2 tt/2 

8\9) 



(82) / W) 4 ^(*)« / 4(W) (sin g )m . 1 rfg+ / 



sin m 3 



34 IGOR WIGMAN 

The contribution of the main term in (|82p to the integral in (|8ip is 

tt/2 

lm-3 



< I 4(N9)-^d9 = N 



o 

l 



'<,<>! . 



km— 3 







ttN 
2 

N m~4 f 

/ d> m - 3 







ttJV 

? # " 




J <\) m - x 


? 


1 




1 1—1 





using the well known decay 

J a (y) < _ 

of the Bessel J functions at infinity. 

The first integral involved in the expression above is O(l), since J a van- 
ishes with multiplicity (at least) a at zero (it follows, for example, from 
Hilb's formula). The second one is bounded by 



log n, m = 2 
1, m > 3 



J./ II l\ I £i 

i 2m ~ A J e m+7 d6 + n~ % J 9 5 ~ m d9 = 0(n m - 12 ). 



Therefore the contribution of the main term in (|82p to the LHS of (|8ip is 
dominated by the RHS of (I5TD . 

The contribution of the error term in (|82p is at most 

l/n w/2 

Jm — 4 

n 

l/n 

We obtain the statement of the lemma by using (f8Tj) and ([73]) with ([78]) . 

□ 

A. 3. Moments of the derivatives of the ultraspherical polynomials. 
Lemma A. 5. 

l 

QTitfd^t) « ^g. 

Proof. We will bound the integral on [0, 1] only, having a similar bound on 
[—1, 0]. By (jT8l) . the statement of the lemma is equivalent to 

l 

J P*'{t) 2 dn(t) «n 2 logn. 
o 

We rewrite the last integral using (|76|) as 

( n + m /2_l)p^ 1 (t)-ntP«(t) 1 



(l_ t 2)2 "^W- 
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To give a bound, we partition the range of the integration into 2 subranges: 



(83) J = J + 

l-l/n 2 

To bound the second integral in ()83|) . we define 

(84) f(t) := (n + m/2 - l)^*) - ntP%(t) = (1 - t 2 )P%'(t). 
Computing the derivative f'(t) and using ([75]) again we obtain 

(85) f'(t) = (m - 2)ti£'(t) - n(n + m - l)i*<(t). 
We claim that this implies 

(86) f{t) « n 2 P«{\) « rT 



»i-2 
2 



the second inequality being a consequence of ([75]) . To see the first inequality 
of (]86p . we note that it is sufficient to show that 



(87) tP«'{t) -^n^. 



by (|85p . With no loss of generality we may assume that t = 1 or i 3 "' has a 
local extremum, i.e. P°"(t) = 0. In both cases the equation ()75[) implies 

i prw = J Pn Q wo(n 2 ), 

which implies (|87|) . 

Now using the linear Taylor approximation of f(t) around t = 1 with 
([56]) ; the second integral in ([83]) is, since /(l) = 0, 

l l 

« n ™+2 J (*~1) 2 2 . (1 _ t2) ^ dt<w m+2 | (1 _^^ dt<<n2 _ 
l-l/n 2 l-l/n 2 



In order to bound the first integral in (j83[) . we employ the generalized 
Hilb's asymptotics (179)) . The integrand is (taking the change of variables 
t = cos 8 and ([ST]) into the account), 



, 2 

1 \ \ DO 



:i + O(£))i£_ 1 (cos0) - cos^(cos« 
n 2 ^ ——3 ^ • (sin0) 2a 

sM.(^ 1 + 0(l))J a ((N-l)e)-[l + 0(9*) + 0(±))J a (N9)j 2 ^ 



■C n ■ — — — ■= h n 



(sin#) 3 (sin#) ; 

2 



J a (N0)-J a {{N-l)0). 2 
« - 2 ~ ¥ L + ^ + O(^0) + n 2 ^, 

and the integration range is essentially f ]. 
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The contribution of the last error term in (|88p is 

vr/2 

< n 2 • rT z / —dO <C 1, 

l/n 



# 3 



the other ones being trivially bounded by 0{n 2 ). 
The contribution of the main term in (188[) is 



T^^-j.aiv-D,)) 2 4 J (ut) - Jam - A))) 



n 2 / ' ~J 3 "-- -^rW * „ ; 



l/n 



<■ n ■ — rd(/) <C n logra, 
n z 7 



2 

1 



due to the boundness of the derivative J' a (t). As it was stated, this is 
equivalent to the statement of the lemma. □ 



Lemma A. 6. One has 

i 



J QT(t) 4 (l-t 2 fdn(t)< 



n 2 log n, m = 2 
^k=z, m>3 



-i 



Proof. The proof of the lemma is similar to the one of lemma IA.51 

We will bound the integral only on [0, 1], having a similar bound on [—1, 0]. 
The statement of the lemma is equivalent to 



(n + m/2 - l)i^Li(t) - nxF£(t) 



l 



n 2 logn, m = 2 



dfi(t) < 



(1 - t 2 ) 2 ^ w I n m , m > 3 

o ^ 

using (f78|) and ([75]) . 

We partition the range of the integration into 2 subranges: 



(89, y = y + 

i_ i 

To bound the second integral in (|89|) we use the definition ([Ml) of the 
function f(t), as well as the inequality ([86]) . as in the course of proof of 
lemma IA.5I Thus the integral is 

l 

« n 2 ^+ 2 ) I - t)^dt « n 2( m+ 2) n - (m+ 4) = n ra_ 

J I 1 "*) 

l-l/n 2 

To bound the first integral in (|89p . we employ the generalized Hilb's 
asymptotics (f79|) . The integrand is (taking into consideration the change 
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of variables t = cos 6), 
(90) 



, 4 

1 \ \ DC* 



;i + O(^))P n a _ 1 (cos0) - co S 0P«(cosl 

4 \ J I ■ /in 4c* 



sin # 3 • sin # 2a 

^.{[l + 0^))U{N-l)e)-{l + 0{6^) + 0{^)J a {N9)j ^ ^ 



■C n • — — — — -: h n 



(sin0) m+1 (sin6») m+1 

4 



Um) -J M N- m , 1 



n 



Qm+l ' 



« n Q^+i + QmTT + 0{n ^=7) + 

and the integration range is (up to a constant) 
The contribution of the last error term in (|90p is 

7I-/2 ?r/2 

< -r / 7; — irdO = tj 2 / r«max n m 4 logra, 1), 

^6 / Qm+l I Qm—1 v ' ' ' 

1/n 1/n 

the other ones being trivially bounded by 0{n m ). 
The contribution of the main term in (|88p is 

(91) 

7T 7riV / -, 

/ (.7.(W)-J.((Ar-iW)' ^ » +4 / ( J ° W - J « W1 - *» 
J o m+l "' J 4> m+1 

1/n 1 

Let g{4>) be the function 

g(t) := J a {4>) - J a (<f>(l - i)). 
Then by the mean value theorem, 

g{4>) = ^ J L(s), 

where 0(1 — jt)) < s < <fi, and using the decay 
we obtain 

(92) \g(d>)\ « ^. 



1 



Substituting (|92|) into (|91j) . we have that the contribution is 



ty 71 7rn 
2 2 



m+4 1 /" 2 m f d(j) _ j n 2 log n, m = 2 



n 4 J (p m+1 J 4> m - x \ n m , m > 3 

1 1 ^ 

which concludes the proof of the lemma. 

□ 
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Lemma A. 7. 

1 

J Q r(t) 2 (l-t 2 fd^t)«-^ 

-l 

Proof. We use the differential equation (j?5l) to write the integral as 
i 

J {mtQ™'(t)-n(n + m-l)Q™(t)) 2 dv(t) 
-l 

l l 
« / {(QT(,t)) 2 d^t)+n 4 J Q r ; i (t) 2 d^t) « -L^ 
-1 -1 
by lemmas I A. II and IA.51 

□ 

Appendix B. The singular functions are "rare" 
In this section we give the proofs of lemmas 12.21 12,31 and 12.41 (see section 

EU). 

Notation B.l. Here and in appendix Owe adapt the following notations. 
Let x and y on the sphere S m such that x ^ ±y. 

(1) Denote xy the (unique) big circle through x and y. 

(2) The smaller arc of xy connecting x to y will be denoted by xy. 

(3) Let z G S m be a point not lying on the plane II = Il(x, y) defined by 
O, x and y. We denote S 2 = S 2 (x,y,z) the (unique) 2-dimensional 
big sphere containing O, x, y and z, i.e. 

S 2 :=S m nU(x,y,z). 

(Note that there is no ambiguity in notations for m = 2). 

We also recall the fact that if S 2 C S m is any big sphere, then for any 
two points x,y G S 2 , 

xy S 2 = xy S m . 

In particular, the shortest path between x and y on S m passes inside S 2 and 

V x d S m(x,y) = V x d s2 (x,y) £ T X (S 2 ) 

under the natural embedding 

T X {S 2 ) C T x {S m ). 

The following simple geometric lemma will prove itself as quite useful. 

Lemma B.2. Let x, y G S m such that x ^ y and £ ^ £' G S m such that 
d(x,£) = d(x,£') and d(y,£) = d(y,£'). Denote v := V x d(x,y) and 

vi = wi(£,0 = V x d(x,£) - V x d(x,?). 

Then for all £ and , v _L v\, and moreover the vectors v\ span v 1 - in 
T x (S m ). 
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Proof. To see the claim of the lemma, we first note that it is obvious for m = 
2. For higher dimensions, it follows from the fact that any 2-dimensional big 
sphere is given by S 2 = S m fl n, where II is a 3-dimensional linear subspace 
of M m+1 , i.e. one direction vector orthogonal to the plane containing xz. □ 

Recall that the set Sing C £ n is the set of singular functions (see definition 

O). 



Proof of lemma \2.2[ We define the map 

f : £ n x S m -► R x R rn 

by 

(/,x)^(/(x),V/(x)), 
using the isometry T x (S m ) = W 71 again, so that 

Sing = vr^(^- 1 (0,0)). 

We claim that \I/ is submersion. Having this claim in our hands would imply 

* _1 (0,0) < M- 1 
by the submersion theorem. Therefore 

dim(Sing) < J\f — 1 

as well. 

To see that ^ is indeed a submersion, we compute its differential to be 

dm = ( m ^ m ^ ••• W ^ *] 

\Vrii(x) Vr] 2 (x) ... Vr/ A r(x) * J ' 

where {%} is the orthonormal basis of £ n , which appears in the definition 
([3]) of /. Denote the matrix j4( m +i) x ./V with the first M columns of . We 
claim that A is of full rank, i.e. rk(A) = m + 1. To see that we compute 
the Gram matrix of its rows to be 



A- A 1 



1 

-zS-Im 



see section [231 Since it is clearly invertible, we conclude that rk(A) = m + 1. 

□ 

Recall that we defined V§ and V% ,y in section O (see (fllj) and (fTSl)). 

Proof of lemma [2~3i . Define B b C SingDV^ to be the set of function having 
±x as their singular point, that is 

Bt = {/ G VI : f(x) = 0, V/(x) = 0}U{/ G £ n : f(-x) = 0, Vf(-x) = 0}. 

It is obvious that B£ is nonempty only if b = 0. Since f(—y) = i/(y) for 
every y G 5 m , / G -B^ implies that / is singular at x. The set i?^ is of 
codimension m, since the covariance matrix (I27p is invertible, so that the 
Gaussian distribution of V/(x) conditioned upon f(x) = 0, is nonsingular. 
Next, we define 

B x b := SingC\V x b \Bt. 
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To prove the statement of the lemma, we need to prove that Bf is of codi- 
mension 1 in Sing n V%. To do so, we follow closely the proof of lemma [2T2l 
This time we define 

^ x :£ n xS m \ {±x} -^R 2 xl ra 

by 

(/,!/)-> (/(*)> /(i/),V/(y)), 

satisfying 

Bf = 7Tfi B (*- 1 (6,0 J 5)). 
Using a similar dimensional approach, it is sufficient to prove that *$> x is 
a submersion. The differential of is 

/ 7}i(x) T] 2 (x) ... w(x) * N 

d^ x = m(y) miy) ■■■ wiy) * 

\Vrfr(y) Vf] 2 {y) ... Vrju(y) *, 

Assume by contradiction that the vectors (r)k(x),r]k(y), ^Vk(y)) satisfy a 
nontrivial linear functional. Since rjk span the whole space £ n , that func- 
tional is satisfied by 

(Q™ (cos d(x, 0), (cos d(y, £)), VQ™(cos d(y, £))) , 

for every £ G S m . The surjectivity of d^> x then follows from lemma |B .31 
Finally, we note that 

which concludes the proof of this lemma. □ 

Lemma B.3. For every x G S m , y ^ ±x, the only solutions in a,/3 G M., 

C G R m for 

(93) 

aQZ(co S d(x,0)+PQ%(cosd(y,0)-Q™Xcosd(y,f;)) S md(y,0(C, V y d(y,0) = 
are a = (3 = 0, C = 0. 

Proof. It is obvious that either a = or (3 = 0, imply that a = (3 = 0, 
C = 0. Therefore we may assume that a = —1, (3 ^ 0. 

Substituting £ in ([93]) and, in addition, any £' 7^ £ not lying on aTy with 
d(£',x) = and d(£',y) = d(£,y), we obtain that C is collinear to any 

Vl = vifa?) G T y (5 m ) of the form 

v 1 = Vyd(y,0-Vyd(y,Z'). 

Lemma IB . 21 implies that C is collinear to v := V y d(x,y). 

We restrict ourselves to any two-dimensional big sphere S 2 C S m contain- 
ing x and y. Knowing that C \\ v, for £ G S 2 on the big circle perpendicular 
to V y d(x,y), ([MD is 

by the spherical cosine theorem, where we denote t := cosd(y,£) and 6 := 
cos d(x,y). It is clear that it is only possible if b = 1, that is x = ±y, which 
is a contradiction. 

□ 
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Proof of lemma \2.4\ To prove the statement of the lemma, we partition the 
set Sing n V b ' y into 

Sing n V^ y = B% U B*' y U 

defining appropriately each of the sets above and proving the statements 
regarding each of them separately. 

First, similarly to the proof of lemma [231 we treat the set 

B x h C Sing n V x h ' v 

of function having ± their singular point, that is 

B x h = {/ G P?" : /(*) = 0, V/(x) = 0}U{/ G £ n : /(-x) = 0, V/(-x) = 0}. 

It is easy to see (exactly as in case of lemma 12 .3|) that B b has codimension 
> 1 in V*' y . 

Next, we treat the case when the function / has a singular point on 
D C S m , a distinguished codimension 1 set of points on the sphere we are 
about to define. Let A\ := xy be the big circle linking x to y and Ai C S m 
be the set of all the points z such that the spherical angle Axzy = § is right 
angle. Define 

D = D x>y = A 1 UA 2 \{±x}. 
It is clear that D is a codimension one set on the sphere satisfying ±x ^ D, 

±y£D - 

Define B b ' y C Sing D V b ' y to be the set of singular functions having a 
-D-point as their singular point, that is 

B = B x b ' y = {/ G V^ y :3zeD: f(z) = 0, Vf(z) = 0}. 

We claim that B has codimension at least 1 in Vf ■ To see that we define 
the map 



by 

It is clear that 



: £ n x D —> R 2 x 



(f,z)^(f(x),f(z),Vf(z)). 



Moreover, is a submersion (see the proof of lemma l2T3j) . 

Therefore, (J& b ' y )~ (pi, 0, 0) is of codimension m + 2 in ^„ x D, i.e. of 
dimension N — 3, so that is of codimension > 1 in Va' y '■ 

Finally, we treat the "generic" case. We define the set 

B = B*' y := Sing n V x h ' y \ {Bf U B**) 

of functions in Sing n 7-^ ,S/ having the set of their singular points outside of 
{±x} U D (i.e. having at least one singular point there). 
We define 

^ y : £ n x S m \(DU {±x}) -^K 3 xR ra 

by 

(f,z)~(f(x),f(y),f(z),Vf(z)). 

It is obvious that 

^ = ^(Kf)- 1 (6i,6 2 ,0,0)). 
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df 



w{v) * 
w( z ) * 



As before, to prove that B is of codimension 1, it is sufficient to prove 
that tyfr' v is a submersion. To see that is a submersion, we compute its 
differential to be 

/ m( x ) mi x ) 
m(y) m(y) 

Its surjectivity follows from lemma IB. 41 
This concludes the proof of this lemma. 

□ 

Lemma B.4. Let x, y and z be points on the sphere S m . Suppose that 
x ^ y, z £ xy, and Zxzy ^ ^. Then the only solution to 

(94) af(x) + 0f(y) + 1 f{z) + (C, V z f{z)) = 

for every f G E n is a = j3 = 7 = 0, C = 0. 

Proof. Substituting 

f(z)=Q™(cosd(z, 0) 

with £ G S m , ([Ml) is 

aQ™(co S d(x, £)) +/3Q™(cosd(y, £)) + Q™(cos d(z, £)) 
-Q™'(cosd(z, £))sind(z, 0(C, V z d(z,0) = 0. 
for £ 7^ ±z. 

Comparing the equality (|95l) for £ not lying on yz and any £' with d(x, £) = 
d(x,r]') and d(z, £) = d(z,£'), we obtain 

Q™' (cos d(z, 0) sin(d(*, 0)(C, V,d(z, £) - V,d(z, £')> 

(96) 



(95) 



P 



g^(cosd(y,O)-Qn(cosd(y,0) 



Let 



with £" 7^ £ be the unique pair of points with 

d(z,e) = d(z,e') = d(z,o 

and 

V z d(z, O - V z d(z, £"') = V z d(z, - V z d(z, O- 
In particular, we have 

d(x,e") = d(x,n- 

Substituting £" and into (f96l) . as we may do, yields 

Q™' (cos d(z, 0) sm(d(z, 0)(C, V z d(z, £") - V z d(z, £"')) 

(97) 







and comparing (|96|) to (|97j) we see that either /? = or 
(98) 

Q™(cos d(y, 0) - <T (cos d(y, ? )) - Q-(cos d(y, £")) + (cos d(y, ?»)) = 0. 
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Suppose by contradiction that the latter holds. We restrict ourselves to 
any big two-dimensional sphere S 2 (x,y,z) C S m (recall notation IB. 1|) . Let 
d > be a small number and (f> ^ (/)' £ S 2 be the (unique) points which 
satisfy 

d(z, 4>) = d(z, (/>') = d 

and (ftp' _L xz. We may approach to <f> by 5 2 -points £ and £" and to <fi' by £' 
and of the form as above with the additional requirement 

d(£,xz) = d(£ ,xz) = d(£ ,xz) = d(£ ,xz) = d. 

Dividing ((98]) by d(£,£") = d(g,g"), and taking the limit as £ -> 0, we 
obtain 
(99) 

(cos d(y, <p)) sin 0) Q^d(y, 0) = Q™(cos d(y, <p')) sin d(?/, <P')-Q^rd(y, <p'), 

where and e^' are the unit tangent vectors in the directions </>£" and / £ /// 
respectively. 

Denote c?o := d(y,z), d\ := d(y,cp) and c?2 := d(y,<f>'). Let <5 be the angle 
5 = Zxzy. We have <5 / 0, ^ by the assumptions of the lemma. We compute 

c\ = c\ (d) := cos d\ = cos do cos d + cos 5 sin do sin c? 

and 

C2 = 02(d) := cos c?2 = cos do cos d — cos (5 sin do sin d, 

by the spherical cosine theorem. It is obvious that the LHS of (I99p is an 
analytic function of c\, and the RHS is the same function evaluated at 
c 2 = g( c i) f° r some analytic function g. The function g is defined on an 
neighbourhood of cos do satisfying g(cosdo) = cos do- Therefore, lemma lB~5l 
implies that 

g' (cos d ) = ±1. 

On the other hand, computing the derivative explicitly, we have 



g' '(cos d ) 



cos do sin do — cos 5 cos do sin do 1 + cos 6 



- cos do sin do + cos 8 cos do sin do 1 — cos 5 ' 

which, clearly, under the assumptions of the lemma, cannot be equal to ±1, 
and therefore we obtain the necessary contradiction. This proves that = 0. 
By the symmetry, we have a = as well. 
Thus ([961) implies that 

C± Vl (0 :=V 2 d(z,0-V z d(z,O, 
for every £, £' of the form above. However, for every £, the vectors 

«i(U')eT 2 (r) 

are all orthogonal to v := V z d(x,z), the vector in the direction of xz, and 
moreover, they span the orthogonal complement v ± , by lemma IB. 21 

Therefore C must be collinear to v. Similarly we may argue that C is 
collinear to v' := V z d(z, y). However, v and v' are not collinear by the 
assumptions of the present lemma, so that C = 0. Knowing that, 7 = is 
easy to obtain. 

□ 
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Lemma B.5. Let f(t) an analytic, not identically vanishing function, and 
g(t) a differentiable function defined on an neighbourhood I of to £ I such 
that g(to) = to. Suppose that we have on I 

f(g(t)) = f(t). 

Then g'(t ) = ±1. 

Proof. We have by the chain rule, 

f'(g(t))g'(t) = /'(*). 

Therefore, if /'(to) / then g' (to) = 1 and we are done. Otherwise, we 
continue differentiating to obtain 

f"(g(t))g' 2 (t) + f'(g(t))g"(t) = f"(t) 

so that if /"(to) 7^ 0, we have g' 2 (to) = 1 and we are done again. Otherwise 
we continue differentiating until we encounter the first derivative f^ k \to) ^ 
implying g' k (t) = 1. Such a number k exists, since / is analytic. 

□ 

Appendix C. Non degeneracy of point value and gradient 

DISTRIBUTION 

In this section we prove that for ±N ^ x G S m , the distribution of the 
random vector Z defined in section 13.11 is nonsingular Gaussian. 

Lemma C.l. Let x ^ ±iV G S m and V = V X be vector space 

V = M 2 x T x {S m ) x T N {S m ). 

Define the subspace 

U = U x , n C V 

by 

U = {{f(x), f(N), V/(x), V/(y)) : / G £ n }. 

Then one has 

U = V, 

provided that n large enough. That is, the distribution of the random vector 

V = (f(x), f(N), Vf(x), Vf(N)) 
is Gaussian nondegenerate and one may identify 

U ^ M 2m+2 , 

as in section \2Ji\ 

Proof. Let x ^ ±N. We assume by contradiction, that U is a proper sub- 
space of V, i.e. there is a nontrivial functional h : V — * R vanishing on 
U. 

We wish to work with coordinates and employ the orthonormal bases for 
T x (S m ) and T^(S m ) chosen in section [27^1 so that under the corresponding 
identification, one has (|26p . 

By our assumption, there exist numbers a, (3 G R and vectors C, D G R 2 
so that 

(100) af(x) + (3f(N) + (C, Vf(x)) + (D, Vf(N)) = 0. 
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We know that for every rj £ S m , the function 

f(x) := Q™ (cos d(x,rj)), 

is a spherical harmonic lying in £ n . For this particular function (|lUUp is for 

77 / ±x,±N, 

aQ™(cosd(x, rj)) + (3Q™(cosd(N, rj)) 

(101) - QT (cos d(x, rj)) sin (d(x, ri)) ■ (C, V x d(x, rj)) 
-QT(cosd(N, ri))sm(d(N, rj)) ■ (D, V N d(N,rj)) = 0. 

First choose rj £ S m not lying on xN and compare (jlOip for 7/ and any 
rj ^ rj satisfying ?/) = d(x, rj) and rj') = d(y, rj). We obtain 

Q™'{cosd(x, ri))sm(d(x, r,)) ■ (C, V x d(x,rj)) 

+ Q™' (cos d(N, rj))sm(d(N, rj)) • (D, V N d(N, V )) 

= QT(cosd(x, rj))sm(d(x, r,)) • (C, V x d(x, V ')) 

+ QT (cos d(N, rj))sm(d(N, rj)) • (D, V N d(N,rj')). 

Equivalently, 

Q™'(cos d{x, rj)) sin (d(x, rj)) ■ (C, V x d(x, rj) - V x d{x, rj')) 

(102) = -Q™'(cos d(N, rj)) sin (d(N, r,)) ■ (D, V N d(N, rj) - V N d(N, r,')). 
For every rj, the vectors 

vi(rj) = V x d(x,r,) - d(x,rj') £ T x (S m ) 

are all orthogonal to v := V x d(x, N), the vector in the direction of xN, 
and moreover, they span the orthogonal complement v 1 - by lemma IB. 2 1 We 
claim that the equality (I102p implies that 

C _L sp{vi(rj)} 

and thus C and v are collinear. Similarly, D and 1/ := V ^d(x,N) are 
collinear, and since we identify v with —v', that implies C = XD are collinear. 

Suppose otherwise. Let vq = vi(i]q) such that (C, vq) / 0, and consider 
the two-dimensional sphere S 2 C S m defined by x~N and vq. For rj £ 5 2 , 
one has 

ui(ry) || f - 

We fix d = d{N,rj) so that cosd is a zero of Q™' . Then the RHS of (fT02j) 
vanishes and our assumptions imply that cos d{N,rj) is a zero of Q™'. 

However the function cosd(x, ry) is a continuous nonconstant function of 
rj on the arc 

A := {rj : (i(iV,7?) = d} C S 2 , 

and therefore its image contains an interval, contradicting the finiteness of 
number of zeros of Q™'. Therefore 

C || v, 

which proves our claim, i.e. C = XD for some A £ 1R. 

Substituting the last equality into (I102p with rj £ S m such that 

d(x,T)) = d(N,rj), 
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implies A = — 1, i.e. 

(103) C = -D. 
Now substitute 77 — ► x in (jlUlj) to obtain 

(104) a + (3Q™ (cos d) + Q™'(oos d) sin (d) • (C, V N d(N, x)) = 0, 
where d = d(x,N). We obtain similarly 

(105) aQ™(cos d) + 0- QT (cos d) sin (d) • (C, V^x, AT)) = 0, 

upon substitution 7] ^ N. Since in our identification, we have V x d(x, N) = 
-V N d(x,N), (P3B|| together with (fl"05D imply 

(106) a = /?, 
since 

(Q^(cos d) ^ 1) <= (d ^ 0, tt) ^ (x ^ ±AT). 

We claim that a = and C = 0. Assume otherwise. Consider any two- 
dimensional sphere S 2 containing x~N, and the big circle E C 5 2 defined 
by 

E = {i| £ S 2 : d(x, 7?) = d(N, rj)}. 

On (fTOTj) is, substituting (fl"03"l) and (fl"06j) 

(107) 

2QQ^(cosd(x,77))+Q^'(cosd(x,?7))sind(x,?7)(C, V N d(N, rf)-V x d(x, rj)) = 0. 

It is clear that the vector 

v = V N d(N, 77) - V x d(x, 77) 

is collinear to Va;d(x, N), which, as we have seen, collinear to C. In partic- 
ular, a = if and only if C = and thus we may assume by contradiction 
a ^ and C / 0. 

Since x 7^ ±Af, the point 77 £ E lying on xN, satisfies 

d(x,v) < I 
and the point rj' = —77 £ E satisfies 

d(x,v) > |- 
Therefore there exists a point r]Q £ E with 

d(x,m) = |- 

Then either Q™(cos d(x, 770)) = or Q™'(cos d(x, 770)) = 0, depending on 
whether n is even or odd. However, the equality (I107p implies then 

Q™(cosd(x,7/o)) = Q™'(cosd(x, m )) = 0. 

This contradicts the fact that Q™ does not have any double zeros, since then 
the differential equation ([75]) satisfied by Q™ would imply Q™ = 0. 

□ 
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